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CHARTS  FOR  THE  CONICAL  PART  OF  THE 
DOWNWASH  FIELD  OF  SWEPT  WINGS 
AT  SUPERSONIC  SPEEDS 
By  Jack  N.  Nielsen  and  Edward  W.  Perkins 


SUMMARY 

Analytical  expressions  have  been  derived  for  the  downwash 
throughout  the  entire  induced  flow  field  of  lifting  triangles  of 
infinite  chord  with  leading  edges  either  in  front  of  or  behind  the 
Mach  cone.  These  expressions  have  been  determined  from  the  line- 
pressure  source  theory  of  R.  T.  Jones  and  the  lifting-triangle 
investigation  of  H.  J.  Stewart.  Based  on  these  analytical  results, 
downwash  charts  have  been  prepared  from  which  the  downwash  field 
may  be  determined  for  any  practical  combination  of  leading-edge  sweep 
angle  and  flight  Mach  number.  These  charts  for  the  lifting  triangle 
of  infinite  chord  are  basic  to  the  solution  for  the  downwash  field 
of  any  finite  swept  wing  wherein  the  contributions  of  the  tips  and 
the  trailing  edges  are  to  be  determined  by  the  method  of  lift  can¬ 
cellation  as  employed  by  Lagerstrom  and  others.  The  conical  part 
of  the  downwash  is  obtained  from  the  charts,  and  the  nonconical 
part  must  be  obtained  by  other  means  such  as  the  lift-cancellation 
method. 


INTRODUCTION 

One  of  the  principal  requirements  for  a  rational  analysis  of 
the  longitudinal  stability  of  aircraft  is  a  knowledge  of  the  down- 
wash  field  behind  lifting  surfaces.  Theoretical  methods  based  on 
lifting— line  theory  for  determining  downwash  for  conventional  lifting 
surfaces  at  subsonic  speeds  are  well  known.  However,  before  satis¬ 
factory  agreement  between  experiment  and  theory  was  obtained,  the 
theoretical  downwash  had  to  be  corrected  for  the  local  effects  of 
the  wake  and  the  vertical  displacement  of  the  trailing  vortex  sheet 
(e.g.,  reference  l).  The  methods  of  reference  1  are  applicable  to 
conventional  airplanes  throughout  the  subcritical  speed  range  by 
use  of  the  Glauert-Prandtl  rule  (reference  2).  No  practical  method 
exists  at  present  for  the  calculation  of  downwash  in  the  supercritical 
speed  range. 
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Several  methods  based,  on  solutions  of  the  linearized  differ¬ 
ential  equation  of  compressible  flow  are  available  for  determining 
the  theoretical  downvash  at  supersonic  speedB.  By  the  use  of  the 
conical  flow  theory  of  Busemann  (reference  3 )t  Lagerstrom  in  refer¬ 
ence  -4  has  developed  analytical  expressions  for  the  downwash  field 
of  certain  uniformly  loaded  lifting  surfaces.  By  superposition  of 
these  lifting  surfaces,  Lagerstrom  has  determined  the  downwash  field 
of  a  flat  rectangular  wing  and  a  flat  trapezoidal  wing.  In  addition 
he  has  indicated  a  method  for  determining  the  effect  on  the  downwash 
field  of  adding  a  trailing  edge  to  a  lifting  triangle  of  infinite 
chord  to  form  a  finite  wing.  Heaslet  and  Lomax  (reference  5)  have 
determined  the  downwash  for  points  along  the  intersection  of  the  . 
chord  plane  with  the  vertical  plane  of  symmetry  behind  a  finite  tri¬ 
angular  wing  with  subsonic  leading  edges.  They  have  also  determined 
an  approximate  solution  for  the  downwash  in  the  vicinity  of  this  line 
at  an  infinite  distance  behind  the  wing.  The  Schlicting  vortex  theory 
can  be  used  to  calculate  the  downwash  induced  by  various  span  load 
distributions.  The  concepts  employed  are  analogous  to  the  Prandtl 
lifting— line  theory  and  as  such  consider  only  the  spanwise  lift 
distribution'  and  hence  neglect  any  possible  effects  of  chordwise  lift 
distribution. 

No  determination  has  been  made  of  the  downwash  throughout  the 
induced  flow  field  of  a  swept  wing  in  supersonic  flow.  The  solution 
for  the  downwash  field  for  a  lifting  triangle  of  infinite  chord  is 
basic  to  the  solution  for  the  downwash  field  of  any  finite  swept 
wing  wherein  the  contribution  of  the  tips  and  the  trailing  edge  are 
to  be  determined  by  the  method  of  lift  cancellation  as  employed  by 
Lagerstrom  in  reference  4.  For  a  finite  swept  wing,  the  conical 
downwash  field  of  the  lifting  triangle  of  infinite  chord  represents 
the  entire  downwash  field  except  for  the  regions  of  influence  of 
the  tips  and  trailing  edge.  Within  these  regions,  it  represents 
an  appreciable  contribution  to  the  downwash. 

The  purpose  of  the  present  investigation  was  to  determine 
analytical  expressions  for  the  downwash  throughout  the  entire  induced 
flow  field  for  flat  lifting  triangles  of  infinite  chord  with  either 
subsonic  or  supersonic  leading  edges.  These  expressions  are  used  in 
the  construction  of  nondimens ional  downwash  charts,  which,  when  used 
in  conjunction  with  the  lift-cancellation  method,  provide  a  practical 
Ta&axis  of  determining  the  downwash  field  for  finite  swept  wings.  Although 
the  analysis  was  carried  out  for  Mq2  =  2,  the  charts  are  presented 
in  a  form  making  them  applicable  to  any  supersonic  Mach  number.  The 
charts  are  subject  to  the  usual  limitations  of  linear  theory . 
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SYMBOLS 

Primary  Symbols 
angle  of  attack 

real  part  of  a  complex  function 
free— 6tream  velocity 


free— stream  Mach  number 

cotangent  of  the  sweep  angle  of  the  wing  leading  edge 

slope  of  the  surface  of  the  wedge  section  in  the 
streanwise  direction 


longitudinal ,  lateral,  and  normal  coordinates  with  the 
origin  at  the  apex  and  the  z  =  0  plane  coincident 
with  the  chord  plane  of  the  wing.  (The  positive 
directions  are  indicated  in  fig.  1-) 

oblique  coordinates ,  x’  =  x-my,  y*  =  y-mx,  z'  ■=  z  //l-®2 

potential  of  the  perturbation  velocities 


perturbation  velocities  in  x— ,  y— ,  and  z— directions, 
respectively 

downwash  function  (w  +  iw) 
harmonic  conjugate  of  w 
complex  variable  (se*®) 
a/(  y/x)2+(z/x)2 


l+^/l— (y/x)2—( z/x)2 


tan- x 


■(z/: 


1 

:)J 


'1-J  1-n? 
1  +  */l— m2 


downwash  single  (-v/V0) 
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E( a/1”®2)  complete  elliptic  integral  of  the  second  kind  with 
modulus  1— s? 

complex  variables 

p(!  j/v/i-s04  )  incomplete  elliptic  integral  of  the  first  kind  with 
modulus  V1-eo4  and  Bine  amplitude  | 


U0  j 

E(l,  incomplete  elliptic  integral  of  the  second  kind  with 

modulus  */1-6o4  and  8ilie  amplitude  £ 

/Vl-(l-Bo4 

U  vi-fi2  J 

C  constant  of  integration 

T  real  part  of  a  complex  variable 

5  imaginary  part  of  a  complex  variable 

p  real  part  of  the  incomplete  elliptic  integral  of  the 

r  first  kind 

F  imaginary  part  of  the  incomplete  elliptic  integral  of 

1  the  first  kind 

E  real  part  of  the  incomplete  elliptic  integral  of  the 

r  second  kind 

E  imaginary  part  of  the  incomplete  elliptic  integral  of 

^  the  second  kind 

k  modulus  of  elliptic  integrals  (Vl-So4) 

k.  comodulus  of  elliptic  integrals  ( V1"*2) 


cn 


Jacobian  sine  amplitude  elliptic  function 
Jacobian  cosine  amplitude  elliptic  function 
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dn  Jacobian  delta  amplitude  elliptic  function 

sn(Fr,k) 

Vcr  sn(Fi,k*) 

\x  incomplete  elliptic  integral  of  first  kind  of  sine 

amplitude  |  and  modulus  k 


Subscripts 

P  value  at  point  P 

Q  value  at  point  Q 

c  value  at  Mach  cone 

vc  value  between  plane  wave  and  Mach  cone 

1  refers  to  complex  variable  £ 

2  refers  to  complex  variable 

r  refers  to  real  part  of  a  complex  variable 

i  refers  to  imaginary  part  of  a  complex  variable 


LIFTING  TRIANGLE  OF  INFINITE  CHORD  WITH 
SUPERSONIC  LEADING  EDGES 

In  the  theory  of  supersonic  conical  flow,  the  downwash  field  of 
a  flat  lifting  triangle  is  symmetrical  about  both  the  chord  plane  and 
the  vertical  plane  of  symmetry.  The  downwash  field  above  or  below  the 
flat  lifting  triangle  at  angle  of  attack  a  is  the  same  as  the  down- 
wash  field  below  a  nonlifting  triangular  wing  of  the  same  plan  form 
with  a  streamwise  wedge-shaped  section  of  half  angle  a,  since  the 
flows  above  and  below  the  wings  are  independent.  (See  fig.  1.)  The 
downwash  field  for  this  nonlifting  triangular  wing  may  be  determined 
from  its  streamwise  perturbation  velocity  field  which  may  be  obtained 
from  the  results  presented  by  R.  T.  Jones  in  reference  6.  The  stream- 
wise  perturbation  velocity  field,  or  the  u-field,  for  the  nonlifting 
triangular  wing  is  the  sum  of  the  u— fields  for  two  line  pressure 
sources  coincident  with  its  leading  edges.  Similarly,  the  downwash 
field  for  the  wing  is  the  sum  of  the  downwash  fields  for  the  two  line 
pressure  sources.  However,  since  the  downwash  field  for  one  line 
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source  is  the  reflection  in  the  xz -plane  of  the  downwash  field  for 
the  other  line  source,  it  is  sufficient  to  determine  an  expression 
for  the  downwash  of  a  single  line  source.  The  complete  downwash. field 
is  then  determined  as  the  sum  of  the  field  given  by  this  expression 
and  the  field  given  by  this  same  expression  with  the  sign  of  the 
lateral  dimension  reversed. 

R.  T,  Jones  in  reference  6  gives  the  following  equation  for  the 
perturbation  velocity  u  at  any  given  point  in  the  field  for  a  line 
pressure  source  in  front  of  the  Mach  cone. 


In  this  equation  |dz/cbc|  is  equal  to  a  as  previously  discussed. 

It  should  also  be  pointed  out  that  the  positive  root  of  the  radical 

Vy'2  +  z'2  should  be  used  in  conjunction  with  the  principal  value 
of  the  inverse  cosine.  In  the  succeeding  analysis,  the  downwash 
field  for  the  upper  surface  of  the  nonlifting  wing  will  be  determined. 
However,  since  the  downwash  field  for  this  wing  is  antisymmetric  about 
the  chord  plane,  it  will  be  necessary  to  reverse  the  sign  o&che 
result  for  application  to  the  upper  surface  of  the  lifting  triangle. 

The  perturbation  potential  cp  at  any  point  P  may  be  evaluated 
as  the  line  integral  of  u  along  the  line  parallel  to  the  x— axis 
through  the  point  P. 


cp  = 


(2) 


for  which  the  lower  limit  xi  corresponds  to  a  point  outside  the 
region  of  influence  of  the  line  source  where  the  values  of  <P  and  u 
are  zero.  The  region  of  influence  of  the  line  source  is  bounded  by 
the  Mach  cone  and  the  plane  waves  from  the  line  source  tangent  to  the 
Mach  cone.  (See  fig.  1.)  Figure  2  shows  two  possible  paths  of  integration 
for  equation  (2).  The  path  corresponding  to  point  P  of  this  figure 
intersects  only  the  Mach  cone;  whereas  the  path  corresponding  to  point 
Q  intersects  both  the  plane  wave  (y,2+ z’2  =  0)  and  the  Mach  cone 
(x2  =  y2+72) .  Consider  first  point  P:  The  perturbation  potential 
is  given  by 

Tp  =  J  udx 

</y2+z2 


(3) 
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The  vertical  velocity  component  is  obtained  by  differentiation 
of  cpp .  Thus 


w=  =  r  a* 

£  bzj - -  dz 

Vy2+z£ 

since  at  the  lower  limit  u  *  0  in  accordance  with  equation  (l). 


In  considering  point  Q,  it  is  convenient  to  evaluate  equation 
(2)  in  several  steps:  first  through  the  plane  wave,  then  from  the 
plane  wave  to  the  Mach  cone,  and  finally  from  the  Mach  cone  to  the 
point  Q.  In  -passing  through  the  tangent  plane,  the  value  of  the 
term  cos-1  ( x ' / f\/j'£+z ,2)  in  equation  (1)  jumps  from  0  to  rt ;  thus, 
the  discontinuity  in  u  is  finite,  and  there  is  no  discontinuity 
in  q>.  Hence,  the  lower  limit  of  the  integrals  may  be  taken  as  the 

,  z  f\J  m2— 1 

value  of  x  at  the  plane  wave,  y/m  +  - — - .  Between  the  plane 

wave  and  the  Mach  cone  the  value  of  u  is  constant  at  —  VQma/ ^y/m2— 1 
so  that  within  this  region  the  potential  is 


Voma  f  y  z«/J 
udx  =  -  — -  x  -  - - — 


7  z 

-  +  — - 

m  m 


Vwc  ■  if  *  Y°a 


Vm2-!' 


This  equation  shows  that  the  vertical  velocity  is  constant  in  the 
whole  region  between  the  plane  wave  and  the  Mach  cone. 

When  the  integration  is  carried  from  the  Mach  cone  to  point  Q, 
it  is  noted  from  equation  (5)  that  the  potential  at  the  Mach  cone  is 


*=  - ‘ 


and  therefore  the  potential  at  point  Q  is 


y2+z2 
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Differentiating  to  obtain  tbe  vertical  velocity 


WQ  =  Y°a  + 7  ^  (9) 

//y2+z5 

Equations  (4)  and  (9)  give  the  value  of  the  vertical  velocity  at 
points  P  and  Q.  Although  both  equations  contain  the  same  integral, 
the  actual  value  of  the  integral  will  depend  on  the  path  of  integration, 
The  integral  is  evaluated  in  Appendix  A  for  each  path;  thus,  for 
point  P 


Vy2+z2 

and  for  point  Q 


8u  .  i  yy'+z2 

dx  =  -0-  cos-1  ■  y- —  1  - 

dz  *  ^y2+z2  Jy •2+z72 


f  —  dx  =  — cos-1  C 

{/y2**2  hz  * 


yy'+z2 
y2+z2  /y/y,2+z 


-  Y0a  (11) 


If  the  values  of  the  integral  given  by  equations  (10)  and  (11)  are 
substituted  into  equations  (4)  and(9),  respectively,  identical  results 
for  the  vertical  velocities  at  points  P  and  Q  are  obtained.  Thus 


wp  =  wQ  = 


COS— 1 


f  Ml  1  —  \ 

\  Vy,2+Z  'T  ' 


Equation  (12)  represents  the  vertical  velocity  due  to  one  line- 
pressure  source.  The  vertical  velocity  due  to  both  pressure  sources 
is  obtained  by  adding  to  this  result  the  contribution  of  the  other 
pressurs  sourc©  which,  is  detenninsd  by"  substituting  y  fov  y  n. 
equation  (12).  In  accordance  with  previous  considerations,  the 
vertical  velocity  for  the  triangular  lifting  wing  is  of  opposite  sign 
to  the  vertical  velocity  due  to  the  two  pressure  sources,  thus 


VQa  f  , 
- - - — j  cos1 

*  l 


_ y(y-mx)+z2 

'y2+z2  ^(y— mx)2  — z2  (m2— 1 ) 


+  COS- 


(y+mxj+z2 


y2+z2 
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LIFTING  TRIANGLE  OF  INFINITE  CHORD  WITH 
SUBSONIC  LEADING  EDGES 

Hie  lift  of  a  flat  triangular  wing  with  leading  edges  swept  behind 
the  Mach  cone  has  been  treated  theoretically  by  H.  J.  Stewart  in  refer¬ 
ence  7 .  In  that  paper  a  differential  equation  is  derived  for  the 
complex  downwash  function,  the  real  part  of  which  is  the  vertical 
velocity.  To  determine  the  vertical  velocity  w,  it  is  necessary  to 
integrate  the  differential  equation  and  determine  the  real  part  of 
the  solution. 

In  the  notation  of  this  paper,  Stewart’s  differential  equation  is 


dw  =  2^q  (l+£)2 _ 

^  mE(  Vrf)  (<;2+s02)3/2  (\s  +  -I 

The  actual  integration  of  this  equation  is  accomplished  in 
Appendix  B,  from  which  it  is  found  that 


(lb) 


W  =  w+iw  = 


2V0as0 


mE(Vl-m2)  (l+s02f 


[  (l+2s02)  F(?,Vl-So4) 

+  E(£,  a/1-s04)  +  F(t),  Vl— sq4)  -  E(tj,  \/L-^")]+  C  (15) 


In  this  equation,  the  complex  variables  £  and  q  are  related  to 
y/x  and  z/x  by  the  following  equations: 


(y/x)  +  i  (z/x) 

1  +«Jl-(y/x)2  -(z/x)2 


(16) 


6=_J_ 

a/£2+3 o2 


(17) 


T)  = 


1 

Vl+So2^2 


(18) 


The  value  of  the  constant  C  may  be  determined  from  the  boundary 
condition  that  w  =  -VQa  on  the  wing. 
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The  vertical  velocity  w  is  the  real  part  of  the  downwash 
function  W  in  equation  (15).  Thus,  it  is  necessary  to  determine 
the  real  part  of  the  elliptic  integrals  appearing  in  the  equation. 

The  methods  used  to  obtain  the  real  part  of  the  incomplete  elliptic 
integrals  of  the  first  and  second  kinds  are  described  in  Appendixes  C 
pnH  d,  respectively.  Using  the  results  of  these  appendixes,  the 
final  expression  for  the  vertical  velocity  is 


w  =  -  27°aS° - -j  (l+2s02)  F(<Ai>  V1-so4) 

mE  1— m2 )  (l+s02)2 


_  , _ T.  (l-^o4)  (TlVUl-^l)  [l-(l-^04)^l3 

+  ^s° }  +  -i-o!  n-a-o^i  ~ 


+  F(V^2>  Vl“s04)  ~  V1_s04) 

(l-sn4)  azfJ\z  ( l-*g )  [l-(l-sn4)^g]  j  (19) 

1— 02  [  l- (  l- Sq4  )^2^ 


On  the  basis  of  conical  flow  considerations,  the  vertical  velocity 
is  constant  along  rays  through  the  apex  of  the  lifting  triangle  and 
thus  depends  only  on  y/x  and  z/x.  In  order  to  determine  the  vertical 
velocity  for  a  ray  defined  by  y/x  and  z/x,  the  real  and  imaginary 
parts  of  I  o-rid  t)  must  be  determined  in  accordance  with  equations 
(16),  (17),  and  (18).  From  tx  and  63.,  the  real  and  imaginary  parts  of 
i,  the  values  of  crx  and  Xx  are  obtained  using  equations  (C12)  and 
(C13).  Similarly,  values  of  a 2  and  X2  a^e  calculated  for  the 
complex  variable  tj  using  these  same  equations.  The  vertical  velocity 
w  for  the  ray  defined  by  y/x  and  z/x  is  then  determined  from 
equation  (19)  using  these  values  of  ffx,  Xx,  cr2,  and  X2. 


RESULTS  AND  DISCUSSION 
Downwash  Charts 

The  equations  presented  in  this  report  permit  determination  of 
the  downwash  at  any  point  in  the  induced  flow  fields  of  lifting 
triangles  of  infinite  chord.  The  equations  have  been  used  to  determine 
downwash  charts  for  a  number  of  such  surfaces.  The  range  of  leading- 
edge  sweep  angle  used  in  the  calculation  is  sufficient  to  include 
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all  practical  configurations,  and  the  number  of  lifting  surfaces  for 
which  the  charts  hare  been  determined  is  sufficient  to  permit  inter¬ 
polation.  The  charts  are  presented  in  figures  3(a)  to  3(h)  and  show 
lines  of  constant  de/da  within  the  Mach  cone  in  a  plane  perpendicular 
to  the  flight  direction.  Since  the  downwash  field  is  conical,  the 
downwash  pattern  is  similar  in  all  planes  perpendicular  to  the  flight 
direction.  Only  one  quadrant  of  the  Mach  cone  is  shown  since  the 
downwash  pattern  is  symmetrical  with  respect  to  both  the  horizontal 
and  vertical  planes  of  symmetry. 

In  determining  the  charts  for  the  various  lifting  triangles, 
the  values  of  de/da  were  calculated  at  a  number  of  points  in  the 
Mach  cone  for  Mq  =  ^2.  The  lines  of  constant  de/da  were  then 
determined  from  cross  plots.  Although  the  calculations  have  been 
made  for  Mn  = in  which  case  de/da  depends  only  on  m, 
y/x,  and  z/x,  it  can  be  shown  by  the  Glauert-Prandtl  rule  that 
de/da  for  any  Mach  number  depends  only  on  m  V^2-!,  (yV^o2- l)/x, 
and  (z-\/Mo  -l)/x.  These  parameters  have  been  utilized  in  the  charts 

of  figure  3  so  that  the  results  are  applicable  at  Mach  numbers  other 
than  y2. 


In  determining  the  value  of  de/da  at  any  point  in  the  field 
for  a  given  finite  swept  wing,  the  value  of  de/da  for  the 
corresponding  lifting  triangle  of  infinite  chord  is  first  determined 
from  figure  3.  If  the  point  in  question  is  not  within  the  region 
influenced  by  the  trailing  edge  or  the  tips,  this  value  will  represent 
the  total  de/da.  However,  if  the  point  lies  within  the  region  of 
influence  of  the  trailing  edge  or  tips,  the  value  determined  from 
figure  3  represents  only  part  of  the  actual  value  of  de/da.  For  a 
finite  lifting  triangle,  the  remaining  part,  which  consists  of  the 
contribution  to  the  downwash  of  the  trailing  edge,  can  be  determined 
by  the  method  of  Lagerstrom  (reference  4).  In  the  more  general  case 
of  a  finite  swept  wing,  it  1b  necessary  to  take  account  of  the  sweep 
angle  of  the  trailing  edge  and  the  effect  of  the  tips  by  a  similar 
method  of  lift  cancellation. 


Variation  of  Downwash  Distribution  With 
Leading— Edge  Sweep  Angle 


tlnn  plots  0f  figure  3  Show  the  changes  In  the  downwash  distribur- 

““A9  8W6ep  aD«le  of  the  lead^  ^ge  is  Increased  with  respect 

2J2 ^Se^h  ,Wh6n^he  leadl1*  ed«°  ving  is  conside^ly 

ahead  of  the  Mach  cone  (large  values  of  de/da  ls  lar^y 

abo?e^rUi  I*®  ISCh  u0nQ  and  the  rate  of  change  of  ds/da  vith  distant 
above  or  below  the  chord  plane  Is  small  for  considerable  distances  £om 
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the  chord  plane.  Between  the  Mach  cone  and  the  plane  waves  from  the 
leading  edges,  which  are  tangent  to  the  Mach  cone,  the  value  of  de/da 
is  constant  and  equal  to  unity .  The  surfaces  of  constant  downwash 
intersect  at  the  lines  of  tangency  of  the  plane  waves  and  the  Mach 
cone.  As  the  sweep  of  the  leading  edge  is  increased  (m,/Mo  -1 
decreases),  the  line  of  tangency  moves  along  the  Mach  cone  toward 
the  chord  plane  and  the  rate  of  change  of  de/da  with  distance  from 
the  chord  plane  increases.  When  m^/Mo2— 1  approaches  1,  the  lines 
of  tangency  and  the  leading  edge  approach  coincidence  on  the  Mach 
cone;  and  when  m^M^-l  is  equal  to  1,  the  surfaces  of  constant 
downwash  intersect  on  this  line  of  coincidence. 

When  the  leading  edge  becomes  subsonio,  there  are  two  qualitative 
changes  in  the  downwash  pattern.  First,  the  surfaces  of  constant 
downwash  intersect  at  the  leading  edges  rather  than  at  the  lines  of 
tangency,  which  no  longer  exist;  second,  a  region  of  upwash  occurs 
between  the  Mach  cone  and  the  leading  edge.  This  region  of  upwash 
increases  in  extent  as  the  sweep  of  the  leading  edge  increases.  In 
addition,  the  rate  of  change  of  de/da  with  distance  above  or  below 
the  lifting  surface  increases  rapidly,  and  the  downwash  near  the 
surface  of  the  Mach  cone  becomes  very  small.  In  the  limit,  as  the 
direction  of  the  3  eading  edge  coincides  with  the  free-stream  direction, 
the  downwash  field  disappears. 


Additional  Factors  Affecting 
the  Downwash 

When  the  charts  of  this  report  are  used  in  conjunction  with  the 
lift-cancellation  method  to  obtain  the  downwash  behind  finite  swept 
wings,  certain  factors  which  have  not  been  discussed  thus  far  must 
be  considered.  In  the  linearized  theory  of  supersonic  flow  the 
assumption  is  made  that  the  trailing  vortex  sheet  is  coincident  with 
the  extended  chord  plane.  This  assumption,  although  valid  for  very 
small  angles  of  attack,  will  probably  cause  an  appreciable  error  at 
large  angles  of  attack  for  which  the  trailing  vortex  sheet  is  dis¬ 
placed  from  the  extended  chord  plane  due  to  the  action  of  the  downwash. 
An  additional  effect,  which  is  important  in  the  immediate  vicinity  of 
the  wake  boundaries,  is  inflow  to  the  wake  resulting  from  the  viscous 
properties  of  the  fluid.  The  influence  on  the  downwash  at  subsonic 
speeds  of  these  two  factors  has  been  investigated  and  methods  for 
predicting  their  influence  have  been  reported  in  reference  1.  However, 
sufficient  experimental  data  are  not  yet  available  to  permit  the 
development  of  similar  methods  for  supersonio  flow. 

Another  factor  that  is  neglected  by  the  linear  theory  is  the 
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curling  up  of  the  -vortex  sheet  behind  the  wine  tins  Thin  arra~+ 
negleehed  in  subsonic  flow  calculations  tol  c^eSionS hi^pect- 
nrSnhT  T*  H°weTer>  for  ^  low-aspect-ratio  wings  which  will 
my  be  Lportant  ^  ^  ^  °f  superaonlc  aircraft,  this  effect 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics. 
Moffett  Field,  Calif. 


APPENDIX  A 

EVALUATION  OF  INTEGRAL  IN  EQUATIONS  (h)  AND  (9) 
The  integral  to  be  evaluated  is 


in  which 


u  =  -  r.p. 


Ypa  m 

*  V®2-! 


COB  1 


C 


X* 

Vy,2+z 


(Al) 


(A2) 


The  value  of  chi/dz  is  determined  from  equation  (A 2)  and  substituted 
into  equation  (Al),  with  the  quantities  x»,  y»,  and  z'  replaced  by 
y-ncxj  and  z  ,  respectively.  Thus 

1 .  Isa  m  r*  - _ _  (A3) 

*  t  (y^ni)s-s2(m£-l)] 

From  the  identity 


MCA  TN  No.  1780 


-  YoaJ  \  m(y2+z2  )~x^+z  V”12-1!- 

l  L  /„  /Z 2J_»1  „  /t2_t2_7.2 


( y  V m2— l*~z  )  V X^~T 


-  tan-1  y2+z2  y~z  1  ~1  \ 

L(y  Vm2~1+Z)  Vx2-y2-z2  ^  Vy2+z2 

Before  substituting  the  limits,  it  is  convenient  to  replace  the 
inverse-tangent  difference  by  the  equivalent  inverse  cosine. 


Vna.  f  , 

I  =  -Q-J  cos-1 


y(y-ms)  +  z2 

yy2+z2  ,/ (y-mx)2— z2(m2— l) 


Jj2+zz 


In  this  equation  the  positive  roots  of  the  radicals  are  to  be  used 
together  with  the  principal  values  of  the  inverse  cosine.  At  the 
lower  limit  where 

x2  =  y2+z2 


(x-my  )2  =  (y-HBx)2-  z2(m2-l) 

the  nrincipal  value  of  the  inverse  cosine  is  zero  for  y  < z/m  and 

I  te  ,5/.'  It  can  b.  she™  fro.  figure  2  that  y<*/«  correspond 

to  point  P  and  y>x/m  corresponds  to  point  Q. 


Therefore,  for  point  P 


2“  dx  = 

ciz 


V  y2+z2  Vy,2+Z*2 


and  for  point  Q 


'y2+z2 


VoaT  ./  yy*  +  z2 

22  dx  =  — cos—1! —  - 

8z  rt  L  K^-^Z+ZZ  Jpz+z  '2 


APPENDIX  B 

INTEGRATION  OF  STEWART’S  DIFFERENTIAL  EQUATION 
In  determining  W  from  Stewart's  differential  equation 
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dW  _  2Y0g _ (1+^2)2 _  (m  s 

« '  -(V^>  (tw)3/Y^+  if 

V  i?/ 

it  is  convenient  to  express  the  downwash  function  as  the  sum  of  three 
integrals 


V  = 


2Y0a  sc' 


mE(  )(l-s04) 


1-s 


f/ 


s02  J  /\/^2+Bo2  V^2so2+1 


+(1-S  n2)2 


(l-spg)2  r  dC 


L/; 


0  '  J  (^802)3/VP^i  302  J  (i£S02+l)°'S  JVW 

(B2) 

Introducing  the  change  of  independent  variable 


I  = 


v£  2+soJ 


the  first  two  integrals  become 

r  dt 


v^2+s02V^2s02+1 


=  F(| ,  Vl-So4  ) 


d£ 


VC? 


U2+s02)3/2  A2s02+1  sQ2  Vl-U-So2)^2 


(B3) 


(B4) 


_ S 

(l— so 


a-j-  F(|,  Vl-So4)  +  -  2M1.  4T  Vl-So4 

So  )  so  ^x-so  > 

(B5) 


Making  use  of  the  transformation" 


Vl+302^2 


(B  6) 


1 6 
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the  third  integral  becomes 


Substituting  these  results  into  equation  (B2)  the  equation  for  the 
downwash  function  becomes 


W  =  w+iw  = 


2V0cc,s0 


mE(  \/l— m2 )  (l+sQ2)2 


(l+2s02)  F(i, V1-so4) 


+  C 


(B8) 


The  constant  C  is  to  be  chosen  so  that  v,  the  real  part  of  W, 
is  zero  on  the  Mach  cone. 


APPENDIX  C 

DETERMINATION  OF  THE  REAL  PART  OF  THE  INCOMPLETE 
ELLIPTIC  INTEGRAL  OF  THE  FIRST  KIND 

The  Jacobian  normal  form  of  the  Incomplete  elliptic  integral 
of  the  first  kind  in  the  complex  plane  is  determined  as  follows: 

^t+16  *t 

F(T+i6,k  )mf  -==-===  (Cl) 

Jo  V 1— k2 

In  determining  the  real  and  imaginary  parts  Fr  and  Fj.  of  thin 
integral,  it  is  convenient  to  introduce  the  Jacobian  sine  amplitude 

elliptic  function. 


sn  (Fr  +  iFi,k)  =  T  +  i8 


(C2) 
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Making  use  of  the  addition  formula  for  the  sine  amplitude  function 

sn(Fr+iFi;k)  sn(*r>k)cn(  iFj,k)dn(  iFj,k)+sn(  iFj,k)cn(Fr,k)dn(Fr,k) 

l-k2sn2(Fr,k)sn2(  iFi,k) 


(C3) 


Substituting  the  following  relationships 

sn 

'  cn( 

1 


(lFiA)  - 

cnCFi^k*) 


cn( iF1,k)  = 


cn(Fi,k* ) 


>  (C4) 


equation  (C3)  becomes 

t 

sn(Fr+iFi,k)  =  sn^r^^)dn(Fi>kt  )+i  sn(Fj,k«  )cn(F1,k«  )cn(Fr,k)dn(Fr,k) 

cn2(Fi,k’  )+k2sn2(Fr,k)sn2(Fi,k' ) 


(C5) 


Separating  equation  (C5)  into  its  real  and  imaginary  parts 

T  =  _ sn(Fr,k)dn(F1,k») _ 

cn2(F1,k’  )+k2sn2(Fr,k)sn2(Fi,kl ) 

&  _  sn(F1,k> )cn(F1^kt )cn(Fr,k)dn(Fr,k) 
cn^F^k*  )+k2sn2(Fr,k)sn2(Fi,k* ) 


(C6) 

(C7) 


Introducing  the  substitutions 


X  =  sn2(Fr,k) 

(C8) 

a  =  sn2(Fi,k») 

(C9) 

and  making  use  of  the  identities 
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cn(Fr,k)  =  \/l— X. 
dn(Fr,k)  =  v/l-k2^ 
cn(Fi,k’ )  =  \/l-0 

dn(Fi,k')  =  v/l— k,2cr 
equations  (C6)  and  (C7)  "become 

t2  =  M  1-fc,gg) 

(l-o+k2!^)2 

52  =  gtl-tf)  (l-X)  (l-k2!) 
(l— a+k?Xcr)2 


(CIO) 


(Cll) 


Simultaneous  solution  of  these  equations  yields 


_  [(14t2+52)  -  y(l-fT2+62)2  -  4t2]  fl+k2(T2+52)  — l+k2(T2+5g)]  2  -  4k2  t2'} 


Ifk^T2 


(r2+82-\) 


a  = 


(t2+62—X)  -  [Xk^T^+d*)  -  1] 


2/_2  .  ~2 


(C12) 

(C13) 


Therefore,  to  find  the  real  and  imaginary  parts  Fr  and  F^,  of 
the  incomplete  elliptic  integral  of  which  the  sine  amplitude  is  r+15, 
it  is  first  necessary  to  determine  X  and  a  from  equations  (C12)  and 
(C13).  Then  from  equations  (C8)  and  (C9)>  it  follows  that 


=  F(y/X,k) 

(Cl4) 

=  F(Vff,k») 

(C15) 
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APPENDIX  D 

DETERMINATION  OF  THE  REAL  PART  OF  THE  INCOMPLETE 

ELLIPTIC  INTEGRAL  OF  THE  SECOND  £IND 

% 

The  Jacobian  normal  form  of  the  incomplete  elliptic  integral 
of  the  second  kind  in  the  complex  plane  is  given  as  follows: 


r+i5  1 . 

,  x  P  i/l-k2i 

E(T+i8,k)  =  /  v_ 


dS 


(Dl) 


In  determining  the  real  and  imaginary  parts  Er  and  E.j_  of 
this  integral,  it  is  convenient  to  introduce  the  transformation 

I  =  sn(  n ,  k)  =  sn^+iii^k)  (D2) 


where  the  real  and  imaginary  parts  of  n  are  hr  and  respec¬ 

tively.  With  this  transformation  equation  (Dl)  becomes 


E(r+iS,k)  =  J 
where  the  upper  limit 


Fr+iFi 


Fr+iFi  = 


t+18 


dn2(n,k)dp 


d£ 


l  Vl-^/v/l-k2!2 

is  evaluated  in  accordance  with  the  methods  of  Appendix  C 


(D3) 


0*) 


In  integrating  equation  (D3)  it  is  convenient  to  perform  the 
integration  in  two  steps.  The  first  integration  is  along  the  real 
axis  to  Fr,  and  the  second  integration  is  along  a  line  parallel  to 
the  imaginary  axis  from  \x  =  Fr  to  =  Fr+lFj_.  Thus  equation  (D3) 
hecomes 


E(  x+i5,k) 


F- 


dn2(ur,k)dnr+i 


Ji 


dn2(Fr+ipi,k)dpi 


(D5) 


/ 
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With  the  substitution  of  equations  (C8)  and  (C9)  of  Appendix  C,  equation  (D9)  becomes 
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E. 


w  l-aCl-k2^) 


(DIO) 


Thus  the  real  part  of  the  incomplete  elliptic  integral  of  the 
second  kind  of  sine  amplitude  t+18  is  obtained  by  evaluating  X. 
and  cr  from  equations  (C12)  and  (C13)  and  substituting  these  values 
into  equation  (DIO). 
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Figure  /. - Wing  and  Mach  wave  system 
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Figure  3. —  Downwash  contours  for  lifting 
triangles  of  infinite  chord. 
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Figure  3.  —  Continued 
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Figure  3. —  Continued 
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Figure  3. - Concluded 


